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CROSS REFERENCE TO RELATED PROVISIONAL APPLICATION 
[0001] This application is based on Provisional Application S.N. 60/237,957 filed 
on October 4, 2000. 



BACKGROUND OF THE INVENTION 
10 [0002] Field of the Invention . The invention is generally related to methods and 

processes for analyzing well production data and for optimizing production of multilayer 
commingled reservoirs and is specifically directed to a methodology for optimizing 
production using commingled performance data and logging information. 

[0003] Discussion of the Prior Art . Field production performance data and multiple 
15 pressure transient tests over a period of time for oil and gas wells in geopressured reservoirs 
have been found to often exhibit marked changes in reservoir effective permeability over the 
producing life of the wells. Similarly, the use of quantitative fractured well diagnostics to 
evaluate the production performance of the hydraulically fractured wells have clearly shown 
that effective fracture half-length and conductivity can be dramatically reduced over the 
20 producing life of the wells. A thorough investigation of this topic may be found in the paper 
presented by Bobby D. Poe, the inventor of the subject application, entitled: "Evaluation of 
Reservoir and Hydraulic Fracture Properties in Geopressure Reservoir," Society of Petroleum 
Engineers, SPE 64732. 

[0004] Some of the earliest references to the fact that subterranean reservoirs do not 
25 always behave as rigid and non-deformable bodies of porous media may be found in the 
groundwater literature, see for example, "Compressibility and Elasticity of Artesian 
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Aquifers," by O. E. Meinzer, Econ. Geol. (1928) 23, 263-271. and "Engineering Hydraulics," 
by C. E. Jacob, John Wiley and Sons, Inc. New York (1950) 321-386. 

[0005] The observations of early experimental and numerical studies of the effects 
of stress-dependent reservoir properties demonstrated that low permeability formations 
5 exhibit a proportionally greater reduction in permeability than high permeability formations. 
The stress-dependence of reservoir permeability and fracture conductivity over the practical 
producing life of low permeability geopressured reservoirs has resulted in the following 
observations: 

[0006] 1. Field evidence of reservoir effective permeability degradation with even 
10 short production time can often be observed in geopressured reservoirs. 

[0007] 2. Quantitative evaluation of the field production performance of hydraulic 
fractures in both normal and geopressured reservoirs have resulted in the observation that the 
fracture conductivity of hydraulically fractured wells commonly decreases with production 
time. 

15 [0008] 3. Multiphase fracture flow has been demonstrated to dramatically reduce 

the effective conductivity of fractures. 

[0009] 4. Pre-fracture estimates of formation effective permeability derived from 
pressure transient test or production analyses are often not representative of the reservoir 
effective permeability exhibited in the post-fracture production performance. 

20 [0010] The analysis of production data of wells to determine productivity has been 

used for almost fifty years in an effort to determine in advance what the response of a well 
will be to production-simulation treatment. A discourse on early techniques may be found in 
the paper presented by R.E. Gladfelter, entitled "Selecting Wells Which Will Respond to 
Production-Simulation Treatment," Drilling and Production Procedures, API (American 

25 Petroleum Institute), Dallas, Texas, 117-129 (1955). The pressure-transient solution of the 
diffusivity equation describing oil and gas flow in the reservoir is commonly used, in which 
the flow rate normalized pressure drops are given by: 

(Pi-Pyvfj/qcand 
30 {P P (Pi)-P P (P W f)}/q g , 

for oil and gas reservoir analyses, respectively, wherein: 
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Pi is the initial reservoir pressure (psia), 

f w /is the sandface flowing pressure (psia) 

q 0 is the oil flow rate, STB/D 

P p is the pseudopressure function, psia 2 /cp, and 
5 q g is the gas flow rate, Mscf/D 

[0011] While analysis of production data using flow rate normalized pressures and 
the pressure transient solutions worked reasonably well during the infinite-acting radial flow 
regime of unfractured wells, boundary flow results have indicated that the production 
normalization follows an exponential trend rather than the logarithmic unit slope exhibited 
10 during the pseudosteady state flow regime of the pressure-transient solution. 

[0012] Throughout most all production history of a well, a terminal pressure is 
imposed on the operating system, whether it is the separator operating pressure, sales line 
pressure, or even atmospheric pressure at the stock tank. In any of these cases, the inner 
boundary condition is a Dirichlet condition (specified terminal pressure). Whether the 
15 terminal pressure inner boundary condition is specified at some point in the surface facilities 
or at the sandface, the inner boundary condition is Dirichlet and the rate-transient solutions 
are typically used. It is also well known that at late production times the inner boundary 
condition at the bottom of the well bore is generally more closely approximated with a 
constant bottomhole flowing pressure rather than a constant rate inner boundary condition. 
20 [0013] An additional problem that arises in the use of pressure-transient solutions as 

the basis for the analysis of production data is the quantity of noise inherent in the data. The 
use of pressure derivative functions to reduce the uniqueness problems associated with 
production data analysis of fractured wells during the early fracture transient behavior even 
further magnifies the effects of noise in the data, commonly requiring smoothing of the 
25 derivatives necessary at the least or making the data uninterpretable at the worst. 

[0014] There have been numerous attempts to develop more meaningful data in an 
effort to maximize the production level of fractured wells. One such example is shown and 
described in U.S. Patent No. 5,960,369 issued to B.H. Samaroo, describing a production 
profile predictor method for a well having more than one completion wherein the process is 
30 applied to each completion provided that the well can produce from any of a plurality of 
zones or in the event of multiple zone production, the production is commingled. 
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[0015] From the foregoing, it can be determined that production of fractured wells 
could be enhanced if production performance could be properly utilized to determine fracture 
efficiency. However, to date no reliable method for generating meaningful data has been 
devised. The examples of the prior art are at best speculative and have produced 
5 unpredictable and inaccurate results. 

SUMMARY OF THE INVENTION 
[0016] The subject invention is an overall petroleum reservoir production 
optimization methodology that permits the identification and remediation of unstimulated, 
under-stimulated, or simply poorly performing reservoir completed intervals in a multilayer 

10 commingled reservoir that can be recompleted using any of various recompletion methods 
(including but not limited to hydraulic fracturing, acidization, re-perforation, or drilling of 
one or more lateral drain holes) to improve the productivity of the well. This invention is an 
excellent reservoir management tool and includes the overall analysis and remediation 
methodology that has been developed for commingled reservoirs. This invention utilizes the 

15 recently developed commingled reservoir system production allocation analysis model and 
procedures described in my copending application, entitled: "Evaluation of Reservoir and 
Hydraulic Fracture Properties in Multilayer Commingled Reservoirs Using Commingled 
Reservoir Production Data and Production Logging Information," Serial No. 09/952,656, 
filed on September 12, 2001, incorporated by reference herein. 

20 [0017] The specialized recompletion techniques that can be used to improve the 

productivity of previously completed individual reservoir intervals in a commingled reservoir 
include but are not limited to coil tubing hydraulic fracturing, conventional fracture and 
matrix acidizing stimulation techniques that use zonal isolation, and re-perforation of the 
individual completed intervals. 

25 [0018] The subject invention is a method of and process for evaluating reservoir 

intrinsic properties, such as reservoir effective permeability, radial flow steady-state skin 
effect, reservoir drainage area, and dual porosity reservoir parameters omega (dimensionless 
fissure to total system storativity) and lambda (matrix to fissure crossflow parameter) of the 
individual unfractured reservoir layers in a multilayer commingled reservoir system using 

30 commingled reservoir production data, such as wellhead flowing pressures, temperatures and 
flow rates and/or cumulatives of the oil, gas, and water phases, and production log 
information (or pressure gauge and spinner survey measurements). The method and process 
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of the invention also permits the evaluation of the hydraulic fracture properties of the 
fractured reservoir layers in the commingled multilayer system, i.e., the effective fracture 
half-length, effective fracture permeability, permeability anisotropy, reservoir drainage area, 
and the dual porosity reservoir parameters omega and lambda. The effects of multiphase and 
non-Darcy fracture flow are also considered in the analysis of fractured reservoir layers. 

[0019] The production performance of horizontal and slanted well completions, 
including both unfractured and hydraulically fractured horizontal and slanted wellbores, can 
be evaluated using the subject invention to also determine the vertical-horizontal permeability 
anisotropy ratio, and effective horizontal wellbore length. Radial composite reservoir models 
can also be used in the analysis procedure to identify the individual completed interval 
properties of a commingled multilayer reservoir with two or more regions of distinctly 
different properties. 

[0020] The flow rates and cumulative production of all three fluids (oil or 
condensate, gas and water) produced from each completed reservoir interval and the 
corresponding midzone wellbore pressure history are obtained using the commingled 
reservoir production allocation analysis model and procedures presented in my 
aforementioned copending application, in addition to the commingled reservoir production 
history record, and production log (or spinner survey and pressure gauge) measurements of 
the well. The identification of water and hydrocarbons can be determined from the 
production log. If the more advanced gas holdup detection and measurement is used in 
combination with the production log, the gas and hydrocarbon liquid production can also be 
determined from the flowing wellstream fluid. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0021] Fig. 1 is an illustration of the systematic and sequential computational 
procedure in accordance with the subject invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0022] The subject invention is directed to a method for optimizing overall 
petroleum reservoir production through the identification and remediation of unstimulated, 
under-stimulated, or simply poorly performing reservoir completed intervals in a multilayer 
commingled reservoir, permitting recompletion using any of various recompletion methods 
(including but not limited to hydraulic fracturing, acidization, re-perforation, or drilling of 
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one or more lateral drain holes). The method of the subject invention provides a reservoir 
management tool and includes the overall analysis and remediation methodology that has 
been developed for commingled reservoirs. This invention utilizes the recently developed 
commingled reservoir system production allocation analysis model and procedures described 
in my copending application, entitled: "Evaluation of Reservoir and Hydraulic Fracture 
Properties in Multilayer Commingled Reservoirs Using Commingled Reservoir Production 
Data and Production Logging Information," Serial No. 09/952,656 filed on September 12, 
2001, incorporated by reference herein. 

[0023] Fig. 1 is an illustration of the systematic and sequential computational 
procedure in accordance with the subject invention. Beginning at the wellhead (10), the 
pressure traverses to the midpoint of each completed interval are computed in a sequential 
manner. The fluid flow rates in each successively deeper segment of the wellbore are 
decreased from the previous wellbore segment by the production from the completed 
intervals above that segment of the wellbore. The mathematical relationships that describe 
the fluid phase flow rates (into or out) of each of the completed intervals in the wellbore are 
given as follows for oil, gas, and water production of the j* completed interval, respectively: 

q oj (t) = qo t (t)foj(t), 
1gj(t) = QgMfgjit), 
q W j(t) = q wt (t)fwj(t), 

where: 

q 0 j is the j* 11 completed interval segment hydrocarbon liquid flow rate, STB/D, 
q ot is the composite system flow rate, STB/D, 

foj is the j* completed interval hydrocarbon liquid flow rate fraction of total well 
hydrocarbon liquid flow rate, fraction, 
qgjis the j* interval gas flow rate, Mscf/D 
j is the index of completed intervals, 

q gt is the composite system total well gas flow rate, Mscf/D, 

fgj is the completed interval gas flow rate fraction of total well gas flow rate, fraction, 

q wj is the j" 1 interval water flow rate, STB/D 

q wt is the composite system total well water flow rate, STB/D 

f W j is the j" 1 completed interval water flow rate fraction of total well water flow rate, 
fraction. 
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[0024] Once the corresponding fluid flow rates in each segment of the wellbore are 
defined mathematically, using the computational procedure of my aforementioned copending 
application, this data is combined with the commingled reservoir production history record, 
and production log (or spinner survey and pressure gauge) measurements of the well to 
5 determine the most effective recompletion strategy. If more advanced gas holdup detection 
and measurement systems are used in combination with the production log, the gas and 
hydrocarbon liquid production can also be determined from the flowing wellstream fluid. 

[0025] Multiple production logs are considered to properly describe the production 
histories of the individual completed intervals in a multilayer commingled reservoir system. 
10 The crossflow between the commingled system reservoir layers in the wellbore may also be 
specified, using the calculation in accordance with the aforementioned application. All 
measured production log information can be used in the analysis, including the measured 
wellbore pressures, temperatures and fluid densities. The pressure measurements in the 
wellbore permit selection of the best-match wellbore pressure traverse correlation to use in 
15 each wellbore segment. The wellbore temperature and fluid density distributions in the 
wellbore can also be directly used in the pressure traverse calculation procedures. 

[0026] The corresponding fluid phase flow rates in each segment of the wellbore 
are also defined mathematically with the relationships as follows for oil, gas and water for the 
n* wellbore pressure traverse segment, respectively. 

20 

n-1 

qon(t) = qot(t) - Z qoj (t) 

j=l 
n>l 

n-1 

qgn(t)=q g t(t)-Z qgj(t) 
j=i 

n>l 
n-1 

qwn(t) =qwt(t) - 1 qwj (t) 
j=l 
n>l 
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[0027] The flow rate and pressure traverse computations are performed in a 
sequential manner for each wellbore segment, starting at the surface or wellhead (10) and 
ending with the deepest completed interval in the wellbore, for both production and injection 
5 scenarios. 

[0028] The fundamental inflow relationships that govern the transient performance 
of a commingled multi-layered reservoir are fully honored in the analysis provided by the 
method of the subject invention. Assuming that accurate production logs are run in a well, 
when a spinner passes a completed interval without a decrease in wellbore flow rate 

10 (comparing wellbore flow rates at the top and bottom of the completed interval, higher or 
equal flow rate at the top than at the bottom), no fluid is entering the interval from the 
wellbore (no loss to the completed interval, i.e., no crossflow). Secondly, once the minimum 
threshold wellbore fluid flow rate is achieved to obtain stable and accurate spinner operation, 
all higher flow rate measurements are also accurate. Lastly, the sum of all of the completed 

15 interval contributions equals the commingled the system production flow rates for both 
production and injection. 

[0029] In the preferred embodiment of the invention, two ASCII input data files are 
used for the analysis. One file is the analysis control file that contains the variable values for 
defining how the analysis is to be performed (which fluid property and pressure traverse 

20 correlations are use, and the wellbore geometry and production log information). The other 
file contains commingled system wellhead flowing pressures and temperatures, and either the 
individual fluid phase flow rates or cumulative production values as a function of production 
time. 

[0030] Upon execution of the analysis two output files are generated. The general 
25 output file contains all of the input data specified for the analysis, the intermediate 
computational results, and the individual completed interval and defined reservoir unit 
production histories. The dump file contains only the tabular output results for the defined 
reservoir units that are ready to be imported elsewhere. 

[0031] The analysis control file contains a large number of analysis control 
30 parameters that the user can use to tailor the production allocation analysis to match most 
commonly encountered wellbore and reservoir conditions. 
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[0032] The composite production log history and the commingled reservoir system 
well production rates or cumulatives are used to compute the individual completed interval 
production rates or cumulatives. The individual fluid phase flow rates can then be 
determined from the specified individual fluid phase cumulative production or vice versa, for 
5 both the commingled reservoir system wellhead production values and also for the individual 
completed interval values. Either the commingled reservoir system well production flow 
rates or cumulative production values may be specified as additional input. 

[0033] Using the fluid flow rates in each wellbore section, the pressure traverse in 
each wellbore segment is evaluated, specifically the wellbore pressure at the top of that 

10 wellbore section, and the temperature and fluid density distributions in that section of the 
wellbore traverse. This analysis is performed sequentially starting at the surface and 
continuing to the deepest completed interval of the well. The fluid flow phase flow rates in 
each wellbore traverse segment are the differences between the commingled system total well 
fluid flow rates and the sum of the flow rates of the individual fluid phases from all of the 

15 completed intervals above that wellbore traverse segment in the well. Therefore the flow 
rates used in the pressure traverse calculations of the topmost traverse segment in the well are 
the total system well flow rates. For the second completed interval, the fluid flow rates used 
in the pressure traverse evaluation are the total system well flow rates minus the flow rates of 
each of the fluid phases in the top completed interval. The wellbore pressures at the top of the 

20 second pressure traverse are therefore equal to the wellbore pressures at the bottom of the 
first pressure traverse. This process is repeated sequentially for all of the deeper completed 
intervals in the wellbore. 

[0034] From this analysis, a complete production history is computed for each 
individual completed reservoir interval. The complete production history data set includes 

25 the mid-zone wellbore pressures and the hydrocarbon liquid (oil or condensate), gas, and 
water flow rates and cumulative production values as a function of production time. This 
also permits the evaluation of user defined reservoir units that consist of one or more 
completed intervals. The reservoir units can be either fracture treatment stages, or simply 
completed intervals that are located close in proximity together, or simply the users 

30 specification of composite reservoir unit production histories. These individual completed 
interval production histories or the composite reservoir unit production histories are then 
evaluated using one or more of several single zone production performance analyses. 
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[0035] Perforation and gravel pack completion pressure loss models may be 
included to directly compute the sandface flowing and shut-in pressures from the wellbore 
and shut-in wellbore pressures for each individual completed interval. Several perforation 
completion loss models are available in the analyses, as well as numerous gravel pack 
5 completion loss models. 

[0036] The quantitative analysis models used herein invert the individual completed 
interval or defined reservoir unit production histories to determine the in situ fracture and 
reservoir properties in a multilayer commingled reservoir system. The results can then be 
used to identify the unstimulated, under-stimulated or simply poorly performing completed 

10 intervals in the wellbore that can be stimulated to improved productivity. Examples include, 
but are not limited to, various forms of fracturing, acidization, or re-perforation. Fracturing 
operations to recomplete the isolated completed intervals requiring production improvement 
can be conducted using conventional fracture stimulation methodology with zonal isolation 
techniques. Examples include, but are not limited to, sand plugs, bridge plugs, packers, and 

15 squeeze techniques, or with the more recently introduced hydraulic fracturing with coil 
tubing. Similarly, acid stimulation of the poorly stimulated completed intervals can be 
performed using conventional acid stimulation methodology and equipment or with coil 
tubing, with zonal isolation when required. Re-perforation of poorly completed intervals can 
also be accomplished by various means including but not limited to wireline and coil tubing 

20 conveyed perforation methods. 

[0037] Economic evaluation of the production enhancement achieved due to the 
recompletion of the underperforming completed intervals of the well can then be performed 
to determine the viability of various possible and practical recompletion techniques. 

[0038] The invention includes the overall reservoir and production optimization 

25 methodology described in my aforementioned application and utilizes every possible piece of 
reservoir, completion, and production performance information available for the well. This 
includes but is not limited to: open and cased hole well log information; wellbore tubular 
goods and configuration; wellbore deviation hole surveys; perforating and gravel pack 
completion information; well stimulation techniques, treatment execution, and evaluation; 

30 production log, spinner survey, and wellbore measurements; surface separation equipment 
and operating conditions; pressure or rate-transient test data; composite system commingled 
reservoir production data; geologic, geophysical, and petrophysical information and 
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techniques for describing the reservoir; periodic reservoir pressure and deliverability tests; 
and the overall well drilling, completion, and production history. The method is extremely 
flexible and permits consideration of all of the existing well drilling, completion and 
production information that is available, as well as any additional data that is newly acquired. 
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